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Abstract 

The ground state energy of the sinh-Gordon model defined on the strip is studied 
using the boundary thermodynamic Bethe ansatz equation. Its ultraviolet (small width 
of the strip) behavior is compared with the one obtained from the boundary Liouville 
reflection amplitude. The results are in perfect agreement in the allowable range of the 
parameters and provide convincing support for both approaches. We also describe how 
the ultraviolet limit of the effective central charge can exceed one in the parameter range 
when the Liouville zero mode forms a bound state. 

1 Introduction 

It is often taken for granted that the short distance asymptotic of a two dimensional rela- 
tivistic field theory is described by a conformal field theory (CFT). This conception leads to 
the working hypothesis that a massive field theory can be considered as a perturbation of its 
limiting CFT (CPT) by a relevant operator (or by a combination of relevant operators) 
The corresponding (typically dimensional) coupling constant determines the mass scale of the 
perturbed model. This simple scheme holds for the most studied perturbed rational CFT's 
and for certain other models like the sine-Gordon or the imaginary coupled Toda field theories. 
The CPT approach also provides a systematic description of the corrections to the ultraviolet 
CFT asymptotics (with certain reservations concerning the non-analyticity in the couplings of 
the vacuum expectation values, see e.g., [2]). This picture is particularly applicable for finite 
size effects, like the Casimir energy [3], where CPT is applied literally and is often convergent 
(see e.g. [HE]). There are, however, field theory models of different type, where the short dis- 
tance asymptotic is considerably more complicated and, up to now, no systematic description 
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in terms of CPT or something similar is known. Conventionally these models can be called the 
"non-compact" ones, since sigma-models with non-compact target spaces are mostly of this type 
and reveal the same (or sometimes more severe) peculiarities of which we're going to discuss 
now. 

The simplest example is the familiar sinh-Gordon model. This model has been studied for a 
long time and is one of the first discovered integrable theories [H [7|. The factorized scattering 
amplitude is one of the simplest possible and the complete set of form-factors of the basic fields 
is known in a very explicit form [HI [9]. Many other characteristics such as the vacuum energy 
and even the vacuum expectation values (one point correlation functions) [10] are known exactly 
(see below for a brief review). In addition, the most general integrable boundary condition has 
a quite simple form and the corresponding factorized boundary scattering admits a complete 
description [IIl[l2]. (This will be recapitulated briefly in section HI). 

However, it was recognized quite a while ago [13] that the short distance asymptotic of this 
apparently simple model is much more involved than the simple CPT scenario pictured above. 
Even if we do not talk about the short distance behavior of the correlation functions, which is 
not yet well understood even on a qualitative footing (see e.g., |9] for some preliminary results), 
the ultraviolet behavior of the Casimir energy behaves in quite a different way from what we're 
used to in CPT. The corrections to the formal c = 1 CFT predictions behave much softer 
than the usual series in appropriate powers of the scale. It was realized [H] that these leading 
soft corrections are mostly controlled by the so-called Liouville reflection amplitude (LRA), a 
quantity of importance in the explicit construction of the Liouville field theory (LFT) [T4l [TS]. 
(For an explicit construction based on the conformal bootstrap see [IE].) Although there are 
serious arguments to believe that LFT plays also an important role in the description of the 
UV asymptotics of other observables, including correlation functions, the finite size settlement 
is probably the one where our current understanding is the best. The relation between the 
UV behavior of finite size energy and the refiection amplitudes in related non-rational CFT's, 
similar to what was first argued in [H] for the sinh-Gordon case, has been observed in other 
integrable 2D models of "exponential interaction", such as SUSY sinh-Gordon [I7|, affine Toda 
systems [18] and the generalized sausage model [T9ll20]. 

In the present publication we report a study of a somewhat different settlement of the 
Casimir problem where, instead of restricting the system to a finite circle with periodic bound- 
ary conditions, we put it to a finite interval with integrable boundary conditions at both sides. 
The ground state energy in this case is measured by means of a modified version of the ther- 
modynamic Bethe ansatz (TBA), the boundary TEA (BTBA) [21]. The UV corrections in this 
case turn out to be related to the boundary Liouville refiection amplitude (BLRA), the bound- 
ary Liouville version of LRA. The system under consideration turns out to be much more rich in 
physics than the periodic circle one, since the boundary conditions provide enough parameters 
to reach physically interesting regimes. But before turning to these interesting topics, let us 
briefiy remind the standard periodic Casimir effect of the sinh-Gordon model to establish the 
convention. 

The bulk sinh-Gordon model is defined by the Lagrangian density 

AinhG = ^ {da4>f + 2/i cosh (260) (1) 

Here is a two-dimensional scalar field, b a dimensionless parameter and /i a dimensional 
coupling constant which determines the scale of the model. In particular the physical mass m 
of the basic (and the only stable) particle A of the model is related to /i as [22j 



7r/x7(6^ 
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where p is another convenient parameter, often used instead of h 

The model is integrable and its factorized scattering theory is completely characterized by the 
AA — i> AA scattering amplitude 

ofn\ sinh^ -isinvrp 

^y^) = ■ UQ < ■ ■ — • 

smh tl + t sm irp 

The knowledge of the scattering theory allows one to apply the TBA to find the finite size 
ground state energy Eo{R) of the model living on a periodic circle of circumference R 

Eo{R) =£R-^ I cosh e log (1 + e-"(^)) dO (5) 

The infinite volume bulk vacuum energy £ is also known exactly [TOl [23] 

2 

8 sm Trp 

and e{Q^ is the solution to the non-linear integral TBA equation 

mR cosh ^ = e + V? * log (l + e"^^^)) (7) 

(* stands for the convolution in Q). The kernel if>{Q) is related to the ShG scattering amplitude 
dH) as 

id 1 4 sin vrp cosh 6^ 

= -2^d^^°^^(^) = 2^cosh2^^-cos2vrp 

In view of Eqs. ([2], (H [6]), the ShG model possesses the weak-strong duality, 6—^1/6 (or 
p — >^ 1 — p), thus the analysis is restricted to < 6^ < 1 (or < p < 1/2). 
It is also convenient to introduce the "effective central charge" Ceff(-R) 

B„m = (9) 

instead of Eq{R). The most important asymptotic part of the effective central charge at i? ^ 
can be described in terms of the "Liouville quantization condition" as 

Ceff = 1 — 24P^ + power-like corrections in R (10) 

where P is the solution of the transcendental equation [T4] . 

Al(P) = vr + 4PQ log(i?/27r) . (11) 

Al(P) is the phase of the LRA 

S'L(P) = -exp(zAL(P)) (12) 

which reads explicitly 

( ..2^^-2^P/fe ^(l + 2^&P)^(l + 2^&-^P) 
S^{P) = - (vr^7(& )) r(l - 2.6P)r(l - 2.6-P) ' ^'^^ 
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Here fj, is the same coupling constant as in Eq. ([I]) and in the Liouville context is called the bulk 
cosmological constant. Explicit arguments leading to the relation in Eq. (fTTll will be given in 
section [2], for the more complicated case of the "open" finite size effects. We mention here only 
that LET can be obtained formally as a kind of "reduction" of the Lagrangian Neglecting 
one of the exponentials in the interaction term, 2/icosh (260) = /xexp (260) + /lexp (—260), we 
are left with the familiar bulk Liouville Lagrangian 

= ^ (9.0)' + /ie^^^ (14) 
The former is known to define a non-rational CFT with central charge 

cl = 1 + eg' (15) 

where Q is yet another convenient parameter 

g = 6-^ + 6 (16) 

and is, for historical reasons, called the "Liouville background charge" . 

In what follows we are going to apply the same idea of pjj to the system on a finite strip of 
length R with appropriate "right" and "left" boundary conditions (henceforth referred to as "1" 
and "2", respectively), and relate the small R asymptotic of the ground state energy £'strip(-R) to 
the "boundary Liouville reflection amplitudes" 5'b(-P|si, S2) in [2l]. At the same time, EsuipiR) 
can be alternatively "measured" through the BTBA. General integrable boundary condition in 
sinh-Gordon model contains two continuous parameters at each edge (see section [4]), so that 
the "open strip" settlement offers, apart from the overall parameter 6, four parameters to play 
with. This makes the problem quite interesting and rich in physical phenomena. We will start 
with some pedagogical reviews on BTBA and BLRA in the first few sections and provide new 
results in later sections. 

The content is organized as follows. In section [21 we describe briefiy the boundary Liouville 
problem and present the explicit expression for the boundary two-point function first given in 
[21], which coincides with BLRA up to notations. Here the singularity structure and the strong- 
weak duality of the theory are manifest in Barnes double-gamma and double-sine function [25] 
whose definitions and useful relations are found in the Appendix. In section [3l semi-classical 
"mini-superspace" calculation [26] is presented, which gives an independent support to the 
BLRA and will feed our intuition in later discussion. 

We begin section [4] with a brief survey of the factorized boundary scattering in the boundary 
sinh-Gordon model [II],[l2|. Section [5] is devoted to the general formulation of the whole four 
parameter "open strip" problem. Here we develop the usual "zero mode dynamics" arguments, 
which relate the UV behavior of Estrip{R) to the "boundary Liouville quantization equation" 
(involving two different BLRAs, S'b(-P|s^, sj) and S'b(-P|S]~, )). In sectionEl straightforward 
form of the related BTBA [21] is presented and its analytic properties are discussed. The 
standard BTBA equation, however, in a certain region of the parameters needs manipulation 
of the singular behavior of the boundary fugacity to improve the slow convergence of numerics. 
It is to be noted in section [7] that BTBA is insensitive to the sign of the boundary scattering 
parameters, whereas BLRA is not. This mismatch is again due to the singular behavior of 
the boundary fugacity in BTBA and appears in other BTBA problems [27l EH [29] as well. 
BTBA is modified by introducing an additional term, relating to the sign change of the "one- 
particle" coupling in the boundary state [TT]. In section [H the BTBA equation is solved 
numerically and is compared to the result of BLRA. Small R asymptotic from BTBA with at 
least one edge "symmetric" (having non-singular fugacity), is found in excellent agreement with 
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the BLRA result when the parameters are away from the singularity domain of BLRA. Inside 
the singularity domain, BRLA shows that Liouville zero mode is not traveling anymore but 
trapped to form a bound state, which in turn make the UV limit of the effective central charge to 
exceed one. BTBA with both edges "asymmetric' (having singular fugacity) supports perfectly 
the BLRA result. Section is devoted to the analytic calculation of the boundary condition 
dependent UV central charge from BTBA and. The result is found to be in agreement with 
the one coming form BLRA providing another confirmation of the conjectured relation between 
the UV and IR parameters. Section [10] is the summary and discussion. 



2 Boundary Liouville reflection amplitude 

Let us consider Liouville field theory ( fT4ll on a strip of width vr, parametrized by the transversal 
("space") coordinate < cr < vr and the "time" t along the strip. The complex coordinates are, 
as usual ^ = a + ir and ^ = a — ir. Conformally invariant right and left boundary conditions 
are described by the action 



r"00 

A= I dr 



(17) 



Mi,2 = Mo cosh(7r6si,2) , Mq = (—^Y^^ (18) 

Vsm vro^/ 



where, as before, the parameter b is related to the LFT central charge (fTSll with (fTGll . Mi and 
M2 are called the "right" and "left" "boundary cosmological constants" as the bulk cosmological 
constant, and are conveniently parametrized in terms of the dimensionless parameters Si and 

Since Mi 2 are real they will be parametrized as follows: When Mi 2 > Mq, Si,2 = '^1,2 with ri_2 
real and positive. When -Mq < M < Mq, Si,2 = i6"Hl/2 + &Vi,2) with -b-'^ '/2 < (Ti,2 < &"V2 
real. When M < — Mq, there are serious reasons to believe that BLFT is not stable anymore 
and we exclude this range from the investigations. 

Let us denote B the space of states of the LFT on the strip. Conformal invariance entails 
the existence of a (single in this case) set of generators L„ which form the Virasoro algebra 

[Lm, Ln] = (m - n)Lm+n + -^i^^ " m)6m+n , (19) 

It acts on B, splitting B into a set of its highest weight representations. The Hamiltonian, the 
translation generator in r, is 

H=-^ + Lo. (20) 

The best way to understand the structure of B is to take the "zero mode" of the Liouville field 

, , dcr , , 

{a)— (21) 



and consider the region in the configuration space where 0o — oo. Both the bulk and 
boundary interaction terms vanish in this region and we are left with the free massless boson 

on the strip with free boundary conditions on both boundaries. It is convenient, therefore, to 
decompose 0((T, r) in the oscillator modes 

(f){a, r) = 00 - 2tVT + ^e""" cos (na) (22) 
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Here 

^ = -'A- (23) 

while the oscillators satisfy 

[am, flnj = — - — (24) 
The Virasoro generators in this "free field" region are combined as follows 

L„ = ^ akttn-k + (2P + inQ)a„ n 7^ 

Lo = 2^a_fcafc + QV4 + 7^' (25) 

fc>0 

It is easy to argue (see [Hj) that the conformal invariance of the boundary theory prescribes 
the wave function of a primary state \l/p of dimension Ap = Q^/i + to have the following 
asymptotic in the region 0o — > — cxd 

"^P = (exp(2P0o) + 5'b(-P|si, S2) exp{—iP(f)o)) |Fock vacuum) . (26) 

It is the coefficient SB{P\sa, Sb) near the "reflected wave" what is called the BLRA. Exactly as 
in the case of "bulk" reflection [Tl|, the boundary reflection amplitude is unitary 



Sb{P\si,S2) = -exp(zAB(P|si,S2)) (27) 

with the phase Ab(-P|si, S2) real at real P. The standard arguments of real analyticity require 
the analytic unitarity 

SBiP\s,,S2)SBi-P\suS2) = l. (28) 



In the usual boundary CFT language the primary state ( 1261) is interpreted as the one created 
by the (juxtaposition if the right and left boundary conditions are different) boundary operator 

B?/,^:''' = exp((Q/2 + zP)0).,., (29) 

Hence, under a proper normalization the Liouville boundary reflection amplitude is directly 
related to the boundary two-point function Db{P\si, S2) = (^B^^g^{0)B^^g_^{l)'^ , introduced and 
found in 



7^/i7(62)62^262^ 



(Q-2/3)/26 



Here Tb{x) = T2{x\b, b^^) and Sb{x) = S2{x\b, b^^) are the standard double-gamma and double- 
sine functions invented by Barnes |25j (see the Appendix for a brief list of definitions and useful 
relations). 

The BLRA is simply the same quantity with /3 = Q/2 + iP 

Sb{P\si, S2) = ^b(| + ^P\su S2) (31) 
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It is a meromorphic function of P and its phase allows the power expansion, 



A^{P\s,,s,) =J2Msi,S2)P"'~' . (32) 

k=l 

For practical calculations it is convenient to represent the phase (i27l ) in the form 

Ab(P|si, S2) = ^Al(P) + A(P|si, S2) (33) 

where Al(P|si, S2) is the bulk Liouville reflection phase (|T2ll . while the s-dependent part admits 
a convenient integral representation 

00 

^ , , f sm(2Pt)dt cos(sit) cos(s2t) - cosh(6t/2) cosh(b-^t/2) cosh(gt/2) , , 

= y 1 sinh(6t)sinh(t/6) ^^'^ 

—00 

where g = b^^ — b. It is essentially a Fourier transform and is very convenient for numerical 
implementation. 

The states with real values of the momentum P generally constitute the continuous spectrum 
of "physical" states 

B = ®p>oBp (35) 

All these states are characterized by the energy E > —1/24. There are, however, certain 
situations when additional "discrete" states appear. This happens if the reflection amplitude 
Sb{P) has a pole at some P = P^ with QmP^ > 0. Then in the asymptotic (1261) . the incident 
wave is absent and the state is localized. The pole can appear when at least one of Mi or M2 
is negative enough. In the a-parametrization the related poles in the BLRA ([3T]) appear at 



Pri = i { — !— — n ] b ; n = 0,1 



, J- , . . . , 



(CTI +0-2-1) 



(36) 



where [a] stands for the greatest integer less than or equal to a and the boundary parameters 
are limited as 

1 < (Tl + (T2 < 6"^ (37) 

(The right hand side inequality follows from the requirement of the stability of the system 
Ma > -Mo). 

For later comparisons let us quote here the semi-classical expression for BLRA ([3Tll : in the 
limit b —>■ 0, P — > and Si 2 00 while k = P/b and o"i 2 kept flxed 



c^C^DrM - /•47r/i\-^ r(2zfc)r(l/2-a-zA:) 

^ ^"^-^ b^ J T{~2^k)T{l/2-a + ^k) ^^^^ 



where a = (ai + cr2)/2. 



3 "Mini-superspace" approximation 

In this section we provide a semi-classical confirmation of the BLRA ( l38l ). Let us consider the 
semi-classical regime where 6-^0 while k = P/b and the boundary parameters ai^2 are kept 
fixed so that 

Mi,2 = -(7r/i)i/26(Ti,2 (39) 
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In the mini-superspace approximation one neglects all the oscillator modes, replacing the Fock 
space by the vacuum state, and takes into account only the dynamics of the "zero mode" (l2Tll . 
The Hamiltonian (l20l) is replaced by 



1 ^2 

^0 



i^ms = - + 7ifie'''^° + (Ml + M,)e'^'' (40) 



The corresponding eigenfunction of momentum P = bk solves the second order linear differential 
equation 

^ + 7r/ie2^<^» + (Ml + M2)e''A V^(0o) = fc'^(0o) • (41) 

on J 

This is a degenerate hyper-geometric equation. Appropriate solution is 

^(0o) = {^7r^,b-y^' ^^^^l^Z^y^^-^ (2(vr/.)V26exp(60o)) (42) 
where a = {ai + /2 and 

II+1/2 -z/2 roo 

Wx,f,{z) = —— -/ e-'H^-^-^/\l + t)^+^-^/^dt (43) 



r(l/2 + /i- A) 
is the Whittaker function [30]. At 0o ^ — oo 



, . . . ^ e^^^^" - ^(l + ^fc)^(l/2 + ^fc)^(l/2-^fc-a) / vr^x -^^ 
^(l-^fc)^(l/2-^p)^(l/2 + ^fc-a) Ufe^; 



Thus the boundary reflection amplitude in this approximation reads 

Q(ci)r ^ = _ f ^ + 'k)r{l/2 + zk)T{l/2 -ik-a) 

UbV r{i-ik)r{i/2-ik)r{i/2 + ik-a) 

in complete agreement with the corresponding limit of the exact BLRA (l38l) . 



(45) 



4 Boundary sinh-Gordon scattering 

In this section we analyze the sinh-Gordon model with the Lagrangian 1^ in the half-space 
y < 0. The boundary theory is specified by the boundary action, which in the most general 
integrable case has the form [11] 



A 



BshG 



y<o L 



■^{da(pf + 2fi cosh(260) 
47r 



d'^x + j [M+e^^{Q,y) + M~e-^^{Q,y)\dy (46) 



It will be convenient to parametrize the boundary coupling constants following ( ITSll through 
the (self-dual) parameters s"*" and s~ as follows 

M^ = Mo cosh(7r5s=^) (47) 

Integrable boundary conditions are characterized either through integrable boundary inter- 
actions or through factorized boundary scatterings. The relevant amplitude of the factorized 
off-boundary scattering A{6)B = R{6)A{—6)B is easily figured out from ref.[II]. It reads as 

m 

i?(^) = i?o(e)i?{i)(^|r^,^) (48) 
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where the "minimal" amplitude Ro{0) is independent of the boundary parameters 



sinh {9/2 + m/A) cosh(^/2 - mp/A) cosh(^/2 p) /4) 

°^ ^~ smh{e/2-t7c/4)cosh{e/2 + i7Tp/4)cosh{e/2 + i7i{l-p)/4) ' ^ ^ 

Note that -Ro(^) is singular at 6* = z7r/2, which corresponds to the emission of a zero momentum 
particle by the boundary state in the crossed channel, see [H] for the details. 
The second multiplier gives the boundary parameter dependence 

^\ sinh 6* — i cosh (p?7) sinh 9 — i cosh.{p{}) 
' sinh 9 + i cosh(pr]) sinh 9 + i cosh(p'd) 

Here r] and are related to the self-dual parameter [24l [31] 

2bri = 7r(s+ + S-) , 2b^ = 7r(s+ - s') (51) 

Henceforth, we will call the "symmetric" boundary the one with M"*" = M~ (or = s~ = s). 
For the symmetric boundary we have 

br] = 7is, 1^ = (52) 

Let us quote here the expression for the boundary energy fijj.'d) as the function of the 
boundary parameters r] and ^ 

Tfi 

f{ri, ■}}) = ^ . ^ — - (2 cosh(p?7) + 2 cosh(pi?) — sin(7rp/2) — cos(7rp/2) — 1) (53) 

where, as in ([61), m is the mass of the fundamental particle of the sinh-Gordon scattering theory. 
At the best knowledge of the authors this expression has never been obtained rigorously. The 
best way to derive it is to apply the standard relation between the BTBA kernel and the bulk 
and boundary energy (see [32] or [3T] for details). However, strictly speaking this relation is 
justified only in the case of a standard ultraviolet pattern of perturbed rational CFT. It is 
the regular perturbative structure of the short distance corrections which allows to require the 
cancellation of the linear and constant terms [M]. In the case of the sinh-Gordon theory this is 
certainly not the case. As we mentioned in the introduction, the ultraviolet structure is more 
complicated and it is not clear for us how to ask for such cancellation against a background of 
much bigger "soft" corrections. Another way would be to relate the exact one-point function 
of the boundary operator exp{b(f))s,s to f{7isb~^,0) [24l[3T]. However this is not a derivation, 
since exactly this relation has been used to figure out the relation ([STll between the Lagrangian 
and parameters of the scattering theory. Although to our conviction there are no doubts about 
that expression ( [53l ) is correct, in the absence of a rigorous derivation the analysis presented 
below can be considered as its important support. 



5 Sinh-Gordon on a strip 

Now we are ready to consider the whole problem of the sinh-Gordon model on a finite strip 
with two different boundary conditions at the right and left boundaries. Let R be the width of 
the strip. Apart from the bulk parameters b and /x, we need four extra boundary parameters 
M^2 to characterize the boundary interaction at the right and left boundaries. They enter the 
strip action 

/oo 
Lstnp{y)dy (54) 
-oo 
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in the following way 

Lsiripiy) = j^ (^^(9,0)2 + 2/icosh(260)jrfx (55) 
+ M+e^^(0) + Mf e-'"^(0) + M+e^'^{R) + M^e-''^{R) . 
Accordingly we need four parameters sf^2 the usual way related to Mf^g 

= Mo cosh(7r6s^2) (56) 

Scaling properties of the bulk and boundary fields allow to reduce the width of the strip i? to tt 
while rendering the R dependence directly to the coupling constants. This is achieved through 
the rescaling x = {R/'7i)a and y = {R/7i)t. The rescaled Lagrangian reads 

r f 1 f R\'^'^'^^' \ 

i^.trip(r) = I -^{da<Pf + 2/i (-\ cosh(260) j da (57) 

+ ( ^ j (M+e'^(0) + Mf e-'"^(0) + M+e'"^{R) + M^e-'"^{R)) 

Notice that the boundary parameters (l56ll are unchanged under this rescaling. 

For our present analysis it is a good idea to single out the "zero mode" (l2T|l of the field 
0(cr, r) and introduce the "oscillator" operators a„ through (l22ll . (l23ll and (l2lll . The following 
Hamiltonian corresponds to (l57l) (compare with (l20l) ) 

— ^strip = - +2 5]^a_fcafc + /i f — j / exp(260(cr,O))rfo- (58) 

- j (M+e'"^(0, 0) + Mf e'"^(0, 0) + M+e^^in, 0) + M2-e-^'^(7r, 0)) . 
Here the exponentials are thought as normal ordered with respect to the operators a„, e.g., 



0) = e^'^^exp I -226^(-l)"^ j exp UibJ^^i-)"] ■ (59) 

\ n>0 / \ n>0 / 



In the present study we are interested in the ground state energy Ei2{R) of the strip system. 
In terms of the Hamiltonian (l58l) it reduces to finding its lowest energy eigenvector ^'o 



if,tHp^O = ^12(i?)^0 (60) 

in the space of states 

-Sstrip = -^2(0o) ® (Fock spacc of oscillators) . (61) 

Of course in general this is a complicated infinite dimensional problem. However, in the narrow 
strip (ultraviolet) asymptotic R ^ 0, which we mostly consider in the present study, there is 

always a wide "free region" if 6|0o| ^ — log (-R/vr)^ max(Mj^2 5 A'-^''^) ) where both the bulk 
and the boundary interaction terms in ( l58ll can be neglected (see |l3] for similar considerations 
in the "closed" geometry). Here an approximation similar to the mini-superspace one of section 
[3] is rightful and 

\I'o ~ (y4+ exp(zP0o) + A_ exp(— zP0o)) |Fock vacuum) 
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where we introduced a (R dependent) "momentum" parameter P. This parameter is fixed by 
the solution of the problem outside the free region, requiring a usual decay of \l/o at (po ±00. 
In our approximation we neglect exp(— 6</)) and exp(— 260) at 0o ^ cxd, reducing the problem 
to the boundary Liouville one. 



TT 



Similar analysis of the interaction at 0o — ^ —00 results in 



5b(P|4,4). (63) 



S^iP\s^,s^). (64) 



Eqs. (I63l) together with ( 1641 1 give the boundary version of the "Liouville quantization condition" 
(BLQC) and for the ground state P is chosen as the solution to the transcendental equation 

- APQ log(i?/7r) + Ab(P|s+ 4) + AB(P|sr, s^) = 2ix . (65) 

Equations (l62l) and (l65l) with the boundary Liouville reflection phases constitute our ap- 
proximation. In the limit i? — > the solution to P is small 

^~-2Qlog(fl/.)- 

Therefore the smaller the value of P, the better is our approximation. On general footings 
we expect that the leading correction to Eqs. fl62l) and (l65l) are of order P?^'^ (see |33] for 
analogous consideration about the cylinder case). In view of (l66l) this means that this correction 
is exponentially small in P and is of the order of exp (— fovr/P). 

6 Boundary TBA equation 

Formally the BTBA equation gives the strip ground state energy [2T] 

E{R) = / cosh 9 log (1 + Ai2(e)e-"(^)) d9 (67) 

4vr J -00 

where e{9) is the solution to the BTBA equation 

6 = 2mR cosh 9 - ip * L{9) ; L{9) = log(l + Xi2e''){9) . (68) 

Here * is ordinary convolution over the real axis of 9. The strip ground state energy is normal- 
ized to compare with the one EuiR) in (l62l) obtained using the BLQC. 

E,2{R) = E{R) +8R + h + h = -^c,s{R) • (69) 

The quantity Ai2(^) is called the boundary fugacity and is given in terms of the boundary 
scattering amplitude in Eqs. (I49l) and (l50l) as: 

Xuie) = K,{-9)K2{9) (70) 
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where -^'1,2(6') = -^^0(^)^1,2(6') is given in terms of the boundary factorized scattering amphtude 

Eqs. (sniEni) 

Ko{e) = Roit7T/2 - 9) , ki^2{e) = i?S(«7r/2 - 9). 
Exphcitly, the fugacity reads 

^ ^^2 ^ ^os'^ 6^ + cos ^ cosh 6^ + sin ^ cosh 9 — cos{r]ip) 

2 cosh 9 — cos ^ cosh 6^ — sin ^ cosh 9 + cos(?7ip) 
cosh 9 — cos{ri2p) cosh 9 — cos('i9ip) cosh 9 — cos(i92 p) 
cosh 6* + cos{ri2p) cosh 6' + cos('i9ip) cosh 9 + cos(7?2P) 

Let us note that A12 is in general singular at 6* = 0, which reflects the one-particle emission at the 
boundary pT|. However, this singularity vanishes if at least one of the boundaries is chosen, 
say, = 0, to be symmetric. Then, the double pole in Kq{—9)Kq{9) is canceled with the 
double zero of ki{—9)ki{9) and the one-particle emission disappears. In this case the numerical 
analysis of the BTBA equation (l68l) does not show any slow convergence. Strictly speaking this 
is the validity range of the original derivation of the BTBA equation (l68l) in [21] and we are 
safe to apply it only in this domain. Its extension for the case when both one-particle boundary 
couplings are non-vanishing requires some care and we devote the next section to this issue. 



7 BTBA equation with both edges asymmetric 

If both boundaries are asymmetric we face with conceptual and numerical problems. As the 
infrared analysis in [3lj showed the BTBA equation (l68ll describes properly the ground state 
energy only for (71(72 > 0. Even in this case the numerical analysis is also in trouble since 
even though the convolution integration in ( i68l l is flnite, numerical evaluation is very slow in 
convergence. 

To avoid this, one may rewrite the BTBA so that the convolution of the singular part is 
analytically integrated out. One way is to put 



where 



e{9) = 2mR cosh 9 - t{9) - ip * 1^(9) (72) 



L0) = log r""^r/_l. (73) 




1 f cosh — cos7r(p + 7) cosh^ + cos7r(p — 7) 

2 \ cosh6' + cos7r(j9 + 7) cosh6' — cos7r(p — 7) 

(cosh 20 — cos27r(p + 7)) (cosh 20 — cos27r(p — 7)) 
(cosh 29 — cos 2txpY 

Here sin 771 = |(7i(72|e^'"^/2 and (71(72 = 2A/lim5i_>o 0^ Ai2(0) is the residue of the fugacity, 
identifled as 



9a = giVa, ^a) = 2\/ cot ^ cot ^'^^^ tan ^ tan ^ , a = 1,2. (75) 
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This form of the boundary TBA results in the energy of the form 

/OO 
dO coshO U{e) . (76) 
-oo 

where the one particle contribution is manifested at large distances. For gig2 > it is in 
accordance with the boundary analog of the Liischer type correction [34] and L^iO) can be 
further expanded in even number of particle contributions, i.e., powers of e"^"*^ for sufficiently 
large volume. 

There are other choices to de-singularize the BTBA. Another useful form [35], [HH [37] is 
given as 

t{e) = 2mR cosh 9 - ((9) - * L,{9) (77) 



where 



L,{9) = log(tanh2^ + tanh2eAi2(^)e-^(^)) 

/ cosh 9 + sin TTp\ f cosh 29 + cos 27rp \ 



^^^^ J /cosny + smvrpx /cosn/y + cos/7rp\ 
V cosh 6* — sinvrp/ V cosh 26' — cos27rp/ 

The energy in this case has the form 

E{R) = \27i+ d9 cosh 9 Ls{9)\ , (78) 

Not only the infrared analysis suggests that the BTBA equation (i68ll cannot be correct for 
any choice of the boundary parameters if ■ -(92 7^ but this can be seen also from the UV 
analysis: The strip ground state energy obtained from BTBA Eq. (l69ll is not sensitive to the 
sign change of since the boundary fugacity in Eq. (TTTl) is not changed. On the other hand, the 
energy Eq. (l62ll from the BLQC is sensitive to the sign change of as we show now. Suppose 
we change "di into —^i, then according to the relation (ISTl) sf turns into and vice versa so 
that the quantization in Eqs. (1651 ) reads 



- 4PQ log(i?/7r) + A^{P\s^, 4) + Ab(P|4, s^) = 27r . (79) 

This change is serious if t^i ■ 7^ (see Eqs. fl33l [3^ ). Thus, raises a serious question which 
one is correct. 

This mismatch is also noted in the context of other different BTBA problems [27l EH [29] . 
It turns out that the source of trouble is the singular behavior of the fugacity and to cure the 
trouble one needs to modify the original BTBA in Eq. ([68ll when "di ■ 0. 

The correct equation can be obtained by analytical continuation in the one-particle bound- 
ary coupling g m a, model-independent way: To initiate, one notes that the double pole of the 
fugacity induces a pair of zero singularity satisfying 

1 + Ai2(^)e-^(') = (80) 

on the imaginary rapidity axis. This can be easily seen at the infrared (IR) limit. In this case 
putting the zero positions 9 = iu and noting e = 2m_Rcosh6', one has for Eq. ([80ll 

Xi2{iu) e-2-«=°^" = -1. 

The double pole structure of the fugacity results in the zeroes at, with a good approximation 

which is exponentially close to the pole at the origin. 
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Figure 1: Integration contour is deformed including the zero singularity positions 



The convolution integral along the real axis is finite. In order to describe the opposite sign 
of g (or 'd) case, one analytically continues the BTBA equation by deforming the integration 
contour and picking up the zero singularity contribution as shown on Fig. [H If one integrates by 
part the convolution term turns the zero of the logarithm argument to a pair of pole singularities. 
Finally one arrives at the compact form of the new BTBA equation, 

e«(^) = 2mi? cosh ^ + log|^^^^-(^ * log fl + Ai2e-^''') {6) (82) 

b[u -\- lU) V / 

where m > is the positive solution of the zero singularity 

Ai2(m) e-^'''(™) = -1, (83) 

The modified energy has the form 

/oo jn 
— cosh^ log (l + \i2{0) e-^''^(^)j . (84) 

One may see the implication of this result easily at the IR limit. The dominant energy becomes, 
with the help of u > in (ISTl ). 



{R) = m sin m - m^^e-"^ + ■ ■ ■ = m^^e""^^ + • • • , (85) 



which fiips the sign of the IR contribution in (1761) . This result is in agreement with the boundary 
analog of the Liischer type correction [34j. Its confirmation in the UV region will be provided 
by comparing its numerical solution with the one obtained from the BLRA in the next section. 



8 Numerical study 

In the previous sections we presented two different expressions for the strip ground state energy; 
one from BLQC Eqs. and the other from BTBA Eqs. (EUEnD or Eqs. ((821 ES El) . These 

expressions are given either as a transcendental equation or as a nonlinear integral equation 
and are not easy to compare using the analytic expression. In this section, we provide the 
numerical study in a variety of parameter range. 

We first note that the BLRA gives the explicit expression of A(P|si, S2) in flSlj) and Al{P) 
in ([T2I [T3I). On the other hand, the ground state energy is obtained through the BLQC Eq. (l65! l 

A(P|s+, 4) + A{P\s^, S2 ) =271 + APQ log{R/n) - Al{P) , (86) 
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which relates P to the scale R. Thus, to compare the two different approaches, we are enough 
to find the relation R{P) using the BTBA. This relation is obtained via the effective central 
charge through Eq. fl69l) . since the corresponding momentum is given as 



Ptba = V(l-Ceff(i?))/24 

once we use Eq. f l62D . Then, the Liouville boundary phase A*^™"^^ from the BTBA is given as 

-Ai(PTBA). (87) 



A™(PTBA|4,4)+A™(^TBAkr,^2) = 2vr+4PTBAQ log 



TBAj 
TT 



Thus, the numerical check is to compare (l87l) with the analytic expression (l34l) . It is noted that 
the bulk expression Ai(P) in (fT2l [T3ll turns out to be in excellent agreement with the TBA 
result even up to the order of P = 1 (see some of the results in Ref. [Tl]). 
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Figure 2: A(P|s, s) v.s. P with real s: Solid lines refer to the Liouville expression, while circles 
represent BTBA result (6^ = 0.8086 is taken). 



8.1 Symmetric case 

We first restrict ourselves to the case with at least one edge being symmetric, i.e., sf = = si 
or = since in this parameter range, we can avoid the singular behavior of the boundary 
fugacity Ai2(6'). The simplest case is when both boundaries are symmetric so that the right 
edge also has $2 = S2 = S2 or = but S2 is not necessarily the same as si. Here we can use 
that brja = nsa {a = 1,2). In this case, Eq. ( l87l ) simplifies to 

A™(PTBAki, S2) = vr + 2PTBAQl0g ' ^A^(Ptba) . (88) 

In Fig. [2] the numerical results for A*^^^"^^ are compared with the analytic boundary Liouville 
expression for the case of two identical boundaries Si = S2 = s at b'^ = 0.8086 with purely real 
values of s > 0. Another plot is presented for imaginary values of s in Fig. [3l 
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p 


A(P,s,s) 




0.06651092551935 


0.41762001693342 


0.41762001693341 


0.07278412235910 


0.45306073688754 


0.45306073688672 


0.08039322566412 


0.49476603900486 


0.49476603897984 


0.08983119563236 


0.54444716822153 


0.54444716749005 


0.10187859377520 


0.60443320709236 


0.60443318659429 


0.11785948248982 


0.67790492787209 


0.67790438878486 


0.14024177942679 


0.76910412419994 


0.76909136428327 


0.17429836730472 


0.88301300229112 


0.88276474709427 


0.23395539043946 


1.02163214862749 


1.01846195305733 


0.36739025335041 


1.15834702793920 


1.14485912927470 


0.74325569097388 


1.19642293140014 


1.21070661865468 



Table 1: Result for A(P|s, s) v.s. P when s = O.SOi and 6^ = 0.8086. 

At P smaller than 0.15 the numerical agreement is impressively good, as it is illustrated in 
Table [H for the example of s = O.SOi. 

The agreement is also quite good up to P ~ 1 where the values of R are already well bigger 
than the correlation length and we expect the power corrections in R to come into play. 
This can be explained by the fact that after the contributions of the boundary and bulk vacuum 
energy are added in Eq. ( l69ll , the power corrections to Ceg begin with a rather high power of R 
(indeed, they are expected to be ~ P^+^* at 6 < 1). 

Different boundaries with si = s and S2 = s', allow to measure the Liouville reflection phase 
A(P|s,s'). In Fig. mthe results are presented for s = 0.5i and 6^ = 0.8086. 
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Figure 4: A{P\s,s') v.s. P: Solid lines refer to the Liouville expression, while circles represent 
BTBA result (s = 0.5i and 6^ = 0.8086). 

8.2 Discrete mode's case 

As far as one of the edges is symmetric, the two approaches, BTBA and BLRA, are in good 
agreement. However, as Si and S2 approach to the critical value Im(si+S2) = Q, where A{P\s, s) 
becomes singular, the agreement fails except at a small region of P, which is seen in Fig. [3] 
when si = S2 = 0.95i and 6^ = 0.8086 (in this case the actual critical value is s = 1.00565i). 
In fact, there is a parameter range where Im(si + S2) exceeds the critical value Q so that the 
BLRA has the pole at imaginary value of P Eq. (l36l) and at the same time, lm{2bria) < nQ 
and lin{2b'da) < nQ for a = 1, 2 so that there is no boundary bound state in the IR boundary 
scattering theory. This range is given by the following conditions satisfied simultaneously: 

Ka^ + ai at + < 1 , a+ + < I ■ (89) 

In this region, the UV limiting value of Cefr exceeds 1 as shown in Fig. [5l Then, a question 
arises: How is Ces in Eq. (I69l) related to the one Eq. ( l62ll from BLQC? 

At first sight, this question seems not to make any sense since the parameters simply violate 
the convergence of A(P|si, S2). As discussed in section[2l this corresponds to the case when the 
primary operator is not reflected at the Liouville potential wall but is trapped inside as a bound 
state, which can be easily given in the semi-classical approximation in section [3l This happens 
when is sufficiently negative, but not too, so still maintaining the stability of the system. 
This suggests that the Hilbert space has the discrete spectrum as well as the continuous one 
(see also [38l [39]). As a consequence, the bound state energy with Hamiltonian Eq. ( l58l) is 
given as 

^ (4, ^i) 
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Figure 5: Ceff v.s. logR: +-marks are given for p = 0.4 (BLRA predicts Ceff(O) = 1.0225) and 
x-marks are given for p = 0.4471 (BLRA predicts Ceff(O) = 1.02729025). 
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Figure 6: log(— e(i?)) is plotted against log(^) for parameters p = 0.4 and p = 0.4471. Solid 
and dashed lines represent results from BLRA, while the marks + and x correspond to the 
BTBA values. 
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Table 2: Ces{R) 





Cnfff-R) 


log(-e(i?)) BTBA 




1 022497070241 340 


-1 3 1 2561 2063335536 




1 022489761509328 


-1 1 874296674954588 




1 022464134631873 


-10 620393772536636 


-24 Q 


1 022373521967475 


-9 359087748545065 


-1 Q Q 


1 022045862147555 


-8 077079274139289 


-1 4 Q 


1 020790366172398 


-6 736968924454419 


-1 2 1 

J. ^ . J. 


1 018739648008074 


-5 923623092706940 


-If) 1 

J. W . J- 


1 015637398421726 


-5 279891692845977 


log(— e(i?)) obtained from BTBA for various values of 




Ccff(-R) 


log(-e(i2)) BTBA 


-39.9 


1.027288737814910 


-13.883502224634348 


-34.9 


1.027284375181948 


-12.526342850041335 


-29.9 


1.027267391791812 


-11.167552023824138 


-24.9 


1.027200890055835 


-9.803579618174346 


-19.9 


1.026935938305069 


-8.423633052788508 


-14.9 


1.025827041718239 


-6.995003046786642 


-12.1 


1.023913306574199 


-6.139867067007655 


-10.1 


1.020921463329385 


-5.473674443519166 



log(^) at J9 = 0.4: 



Table 3: Ces{R) and log(-e(i?)) obtained from BTBA for various values of log(^) at p = 0.4471 



where P„± is given in ([36] 



n± 



n±]b, n± = 0, 1, 



(90) 



This will give the UV limiting value of the effective central charges greater than 1. 

<±(0) = 1 + 24|P„±|2>1. (91) 

The Ceff(-R) is given in Table [2] and Table [3] and is plotted in figure [5] corresponding to the 
parameters = = 17.5/40,0"^ = 25.5/40, cr2^ = 9.5/40, which satisfies the bound in 
Eq. ((Hi). Eq. (EID predicts c+ (0) = 1.0225 when p = 0.4, and c+j(0) = 1.02729025- ■ ■ when 
p = 0.4471, which agree with the BTBA results. Since BTBA is derived from the saddle point 
of the partition function it always reproduces the lowest energy state. 

We checked this UV limiting value of the effective central charge from BTBA for various 
ranges of the parameters and found a complete agreement. Interestingly, the leading corrections 
to the effective central charge CefT(-R) are no longer logarithmic in the volume but powerlike. 
They are not of perturbative origin, however, but are governed by the analytical continuation 
of the BLQC as we now show. For this we rewrite the BLQC in the exponentiated form 

SB{P\st, st)SB{P\s^, s^) = i-j (92) 

The discrete mode corresponds to the pole singularity of one of the BLRA say S'b(P|s5^, s^). 
As the volume decreases P gets close to the pole at P^+as P = P„+ + ie, so we approximate 
the BLRA in the neighborhood as 

5^(pi4,4) = ^^^^ + ■ ■ ■ = ^^^ii^ + . . . 

P — Pn4- 6 
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For small enough R we determine e from (l92l ) 



/ D\ 4|P„+|Q 

e{R) = G{st, st)SB{Pn+\si, s^) {-] (93) 

The corresponding central charge can be written as Cefr(-R) = l + 24|Pp = 1+24^|P„+| +e(i?) j , 
which gives 

<R) = \/^^^-\Pn,\ (94) 

In the same spirit we compared the BTBA with the exact BLRA we can calculate e{R) from 
BTBA and compare with the expression (1931 ). The data of log(— e(i?)) from BTBA for various 
values of log(^) at p = 0.4 and p = 0.4471 is given in Table [2] and Table [3] and the log-plot 
is given in Figure [8^ The expected slope is 4|P„_|_|Q which is 0.25 for p = 0.4 (and 0.27129 
for p = 0.4471). The fitted value using the lower 4 points is 0.2510 for p = 0.4 (and 0.2719 for 
p = 0.4471). Small R results give more accurate slope and one can see the complete agreement. 
It proves the correctness of the BLRA not only for real but also for imaginary values of P. 



8.3 Asymmetric case 

Next, we are considering the case when both of the edges are asymmetric. The BLQC is given 
as a combination of different boundary Liouville phases. Lets us suppose that -i^i ■ > 0, thus 
we have 

A(P|s+, 4) + A(P|s^, S2 ) =271 + APQ log(P/7r) - Az.(P) , 
then if we switch one of the sign of "i^'s so that {}i ■ '1^2 < 0, we have 

A(P|s+, s^) + A(P|sr, 4) = 2vr + 4Pg log(P/7r) - Ai(P) . 

Since individual phases are confirmed already using the case with at least one boundary sym- 
metric, this settlement looks not to provide any new information to the Liouville phase. Nev- 
ertheless, this combination is important to check the correctness of the analytically continued 
BTBA. 

The numerical results are in perfect agreement with the exact Liouville amplitude if the 
improved version of BTBA Eqs. ((721 [761) or Eqs. ((771 EH]) is applied when ■ ^2 > 0, and the 
modified BTBA Eqs. ((82l [83l El) is applied when 1^1 ■ ^2 < 0. This result is plotted in Fig. [7] 
when Si = $2 = 0, for simplicity. In this case, the phase combination becomes 

^ ^ A(P|.r.0) + A(m4) _ ^ log (^) - Mf) fo, 

z y vr y / 



2 

Other convincing numerical checks are presented in Table ((H) and Table 
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Figure 7: Plot of As{P) (upper black line and green dots, for ■ ^2 > 0) and A^(P) (lower 
red line and blue dots, for -di ■ > 0) for the asymmetric boundaries with = 0.8086 and 
st = S2 = i2b/3 and Si = S2 = 0: Solid lines refer to the Liouville expression and circles to 
the BTBA and modified BTBA result. 



p 


AsiP) 


A^'^A(P) 


0.077653999816446 


-0.039568923372291 


-0.039568923376408 


0.086205185262923 


-0.043897517771010 


-0.043897517774418 


0.096861805449936 


-0.049277259792747 


-0.049277259789521 


0.110502372967937 


-0.056134917955759 


-0.056134917808053 


0.128563844490344 


-0.065154112175326 


-0.065154108313167 


0.153553843986511 


-0.077485060440331 


-0.077484967209008 


0.190215910849849 


-0.095154635137037 


-0.095152687636201 


0.248446714967505 


-0.121766433889740 


-0.121735638152979 


0.351196987049572 


-0.162583158000073 


-0.162317800581419 


0.556084052286610 


-0.216307776233087 


-0.216228485970792 


1.013529248959518 


-0.252173074348906 


-0.257965815998832 


lit for As{P) v.s. P when 6^ = 0.8086 and s+ = s+ = i2b/3 and t 


P 


Aa(P) 


a;^-(p) 


0.078277417261084 


-0.065157660050461 


-0.065157660061779 


0.086965071942987 


-0.072051773672923 


-0.072051773682382 


0.097806408331036 


-0.080511089260293 


-0.080511089203979 


0.111704015386026 


-0.091098590825966 


-0.091098589162912 


0.130133954494199 


-0.104648454630161 


-0.104648417456781 


0.155666785945641 


-0.122400071889728 


-0.122399351214232 


0.193135842394272 


-0.146101177532259 


-0.146089809923920 


0.252485260446018 


-0.177599477043236 


-0.177471319390988 


0.356217516450171 


-0.215816508521331 


-0.215038322905223 


0.560049812334757 


-0.247974299790094 


-0.247287310445203 


1.014342456385944 


-0.256659409495335 


-0.262840890495731 



0. 



Table 5: Result for A^(P) v.s. P when 6^ = 0.8O86 and sf = s+ = i2b/3 and = = 0. 



21 



9 Calculation of the UV central charge from BTBA 



In this section we analyze the UV behavior of the BTBA equation (l68ll . We are able to describe 
the leading small volume behavior of the central charge analytically using the idea developed 
in [I3]: We expand the Fourier transform of the BTBA kernel 

^= /e-VW^ = ^^^^i=i^ (95) 
2tx J ^27r cosh(f ) ' (2n)! 2 ^ ' 



where a = T{^Jp{l — p) = ^ and write the BTBA equation fl68ll for the rescaled functions 
9 ^ 9 + X with X = log mR in the form of an infinite order ordinary differential equation as 



e 



+ e'^-' + log(l - e-^(^)) = log Xu{9) - g (i^^^'"^^) (96) 



We approximate this equation in various rapidity domains in different manners. Since L{9) is 
even we restrict the considerations to the 9 < region. Furthermore we are interested in the 
i? — > , X —>■ —oo limit that is we neglect the e^^~^ term, keeping in mind that we have the 
same contribution from this term in the 9 > domain. We distinguish three rapidity regions 
as follows: x^9,x<^9<ti0 and 6' ~ 0. If 6^ ^ a; the fugacity term is not relevant but we have 
to keep all the derivatives of L. If x ^ ^ we can additionally neglect the term together 
with higher derivatives of L. In this domain, which we call the plateaux domain, L is large and 
positive so we can approximate the BTBA equation as 



yL"(^) + e-^(^) =0 (97) 



The corresponding solution is 



m = io, zr' (98) 



sin^ X{9 

with two arbitrary parameters, A, a, which can be fixed from the boundary conditions. The 
corresponding central charge is 



Ceff (a;) = 1 — + ... (99) 

77^ 



As we decrease the volume, i? ^ 0, the plateaux region becomes larger and larger and the 
approximation is better and better. So in this way we describe the leading correction to the 
central charge. 

The boundary condition at x is provided by the kinetic term as L{x) = 0. In contrast, 
the boundary condition at the origin is determined by the boundary fugacity. If 

/oo 
d9\og{Xi2)<0 (100) 
'OO 

then we can demand the L(0) = — oo boundary condition. This results in 

A = 7 r ; a = 

[a — x) 

and gives the leading UV behavior of the central charge 

(x — a)2 Q^x^ 
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which is in accord with the result of BLRA ( 1661 ). 

In the opposite case when {\ogXi2){k = 0) > 0, the parameter A turns out to be imaginary 
X = in and we have to fit the parameters of the following function: 

, sinh^ k(6' — a) 

m = ^og — — 101 

Demanding L{x) = we have a = x — a. The variable k is determined from the boundary 
condition at the origin, which is provided by the boundary fugacity. Clearly L behaves as 
L oc — 2^6* at the middle of the plateaux region, and we will determine the value of k by 
comparing to the solution around the origin 9^0. Here the kinetic term can be neglected (but 
not the fugacity) and we arrive at the equation 

e = -ip*L , L{e) = log {1 + Xu{0)e-'^^^) (102) 

If we additionally suppose that e is large negative, which follows from L oc —2k9 we arrive at 
the equation 

log Ai2 — L = —(f * L 
which can be solved by Fourier transformation 

(hgXum 



m 



27r [sinh ^ + sinh ^(1 - p) + sinh sinh ^{l-^)-2 sinh ^(1 - ^ 



k 2 sinh ^(1 -p) sinh 

where we used the formula 

dk _.fcg27rsinh^(2-x) 



log([a;]^_e W-+e) = J 



k cosh 



A; 71 



valid for 1 < X < 3 and can be extended for x < 1 via the relation [xjnr g [x]hl+6i 
[2 — x]i2i_g [2 — x]i^^g. We also put rji = ri2 = r] and = ^2 = ^- If in the Fourier transform 

we have a singular term around the origin as — -p then in its inverse Fourier transform we have 
a behavior as —^0. So by inspecting the singularity structure around the origin we can extract 
that 

which no longer depends on x and increases the central charge. The central charge calculated 
from K yields 

Ceff X = 1 + —T- + ■ ■ ■ = 1 + 6^ , + ■ ■ ■ > 1 104 

71'^ P[l — p) 

which agrees with the result coming form the BLRA (iOTI l when the Liouville zero mode is 
trapped in the Liouville potential. This gives a convincing analytical support for the UV-IR 
relation (ISTll and shows the correctness of both the BLRA and the BTBA. 
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10 Discussion 



We have analyzed the ground-state energy of the sinh-Gordon theory defined on the strip subject 
to integrable boundary conditions in two complementary ways using BTBA and BLRA. 

BTBA, being a nonlinear integral equation, systematically sums up the finite size corrections 
to the infinite volume ground-state energy by taking into account the information on the semi- 
infinite boundary scattering theory. As a consequence it is formulated in terms of the boundary 
refiection factors and is reliable in the IR regime. In the case of the boundary scattering theories 
corresponding to perturbed rational BCFTs the careful analysis of the UV limit of BTBA allows 
the determination of the central charge together with the perturbative power-like corrections. 
We have shown in the paper that, in contrast to this usual behaviour, the UV limit of the 
boundary sinh-Gordon theory is governed by a non-rational BCFT: the BLFT. The ground- 
state energy acquires soft (logarithmic) corrections in the volume determined by the BLRA, the 
most important quantity in the bootstrap solution of the BLFT. This approach is valid in the 
UV regime and describes the ground-state energy in terms of the parameters of the Lagrangian. 

As a first step we solved numerically the BTBA and compared with the predictions com- 
ing form the BLRA. In general, we found a convincing evidence of the correctness of both 
approaches. In particular, we checked the previously conjectured relationship between the IR 
and UV parameters (ISTl) and confirmed the predictions of BLRA. Then we used the results of 
BLRA to check the analytically continued BTBA. 

The semi-classical picture, provided in the paper, suggested the existence of a discrete part 
of the Hilbert space, which corresponds to the case when the Liouville zero mode is trapped in 
the boundary Liouville potential. We confirmed the adequacy of this picture at the quantum 
level by numerically calculating the effective central charge, which exceeds one in this case. By 
adopting a method to compute analytically the leading behaviour of the UV central charge we 
were able to derive its value exactly. This provides another support for both the UV-IR relation 
and BLRA. 

Besides confirming the validity of BLRA, which is a widely used quantity in 2D quantum 
gravity, we provided evidence for the discrete part of the Hilbert space. It would be interesting 
to analyze further its consequences. 

The way we performed the analytical continuation in the one-particle boundary coupling 
in BTBA makes possible to apply the result directly to other theories, like boundary Toda 
theories, where the integrable boundary conditions form a discrete set and there is no room 
for playing with any continuous parameter. As a consequence, the refiection factors computed 
from the boundary bootstrap principle in the IR can be compared via the modified BTBA 
to the parameters of the ATFT valid in the UV. This will help to find the sofar unrevealed 
correspondence between the two sets of integrable boundary conditions. 

In |40) the finite volume description of the sinh-Gordon model originating from an integrable 
lattice realization was analyzed. It would be nice to perform a similar analysis for the boundary 
sinh-Gordon theory and explore the analogue of the trapped Liouville mode on the lattice. 

Finally, we note that similarly to the analysis of the bulk staircase model [13] its boundary 
version can be investigated further, along the line of [H], in order to understand better its UV 
limiting theory. 
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Note added This work was almost finished last year while two of the authors visited 
Montpellier and appeared in a final form a month ago. After a communication through the 
phone with Alyosha on October 15th, we were waiting for his last comment on section 9 which 
would be done by Thursday (18th). But very sadly and unexpectedly, we are told he passed 
away on that night. Remembering how excited Alyosha was about the result of this non- 
compact field theory, we put this draft on the web. We have enjoyed life and discussions with 
Alyosha. 



11 Double gamma and double sine 

The double gamma function Ti,{x) was introduced by Barnes [25] through the analytic contin- 
uation in z of the double zeta-series 

oo 

log(b{x,z)= ^ [x + mb + nb^^) ^ (105) 

m,n=0 

convergent if z > 2 (we suppose that 3?e6 > 0). The analytic continuation can be achieved by 
the following integral representation 

where the contour C goes from +oo to +oo encircling the brunch cut of (— t)^ counterclockwise. 
The double gamma function is defined as 



d 



(107) 



z=0 



Like ordinary gamma function, Ti,{x) is a meromorphic function with no zeros and simple poles 
located at x = —mb — nb~^ with (m, n) a pair of non-negative integers. All these poles are 
inside the "wedge" 

|argx| > TT — aigb 

(we imply here that ^mb > 0), which for real b shrinks to the negative part of the real axis. 
Outside the "wedge" it can be represented as the integral which follows directly from (|106p 

f {Q-2xr b^ + b^\ 1 f e-^-log(-t) dt 

where Ce = — r'(l) is the Euler's constant. 

The following dual shift relations are readily derived e.g. from the integral representation 

r^(^ + ^) = y/L;r(6a:) ^^^^^ (109) 
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At large outside the wedge the Stirling asymptotic expansion applies 

(fc-l)!4+2(Q) 



\ogTb{x) 



Mogx + ^ - + 2^ 



24 2 
Here dk{Q) are polynomials in Q defined as 



k=l 



(110) 



4(g) T^nBk-n ^2n-k 



n=0 



n\{k — n)\ 



(111) 



and Bn are usual Bernoulli numbers. One of the eff'ective numerical algorithms is to use several 
times one of the shift relations (whichever is more convenient) to render the argument to a 
region where the Stirling formula with a reasonable number of asymptotic terms is effective. 
For moderate values of x the following form also gives quite accurate results 



logr;,(x) 



g2 _ 2 (g/2 - x)2 \ Ce (g - 2x)2 h'' + b-' 

log X + — 



24 
1 

TT 



4 



12 



^tiWzf log((l + z2)2) 



dz 



(1 - e-^(i+-)V-)(l - e-b-Hw^)Vx^ (i _ iz) 
There is also a convenient line integral representation 

oo 

(X - g/2)2 62 ^ ^-2 



loe Th(x) 



dt 
T 



-xt 



1 g/2 - X 

(l-e-''*)(l-e'*/'') ~^ 



t 
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(112) 



The diperiodic sine ^^(x) (aka as the Barnes double sine function) is related to Tb{x) as 

^'^^^ - UQ - X) 
In the strip < 9fJe x < g it allows the following integral: 

oo 

sinh(g-2x)t (g/2 - x) 



logS'6(x) 



dt 
T 



2sinh(6t) sinh(t/6) 



t 



(113) 



Being a Fourier transform this representation is convenient for numerical calculations. Outside 
the strip of convergence ^^(x) is restored via one of two dual shift relations (which probably 
inspired the name of the function) 



Sb{x + b) = 2sm{nbx)Sb{x) 
Sb{x + 1/b) = 2sm{nx/b)Sb{x) 



(114) 



It is a meromorphic function of x with poles at x = —mb — nb with m and n non-negative 
integers. The only zeros at x = g + mb + nb~^ are predicted by the "unitarity relation" 



Sb{x)Sb{Q - x) 



(115) 



which is a direct consequence of (11131) . The following argument doubling relation is useful to 
arrive at (130]) in the main body of the paper 



Sbi2x) = Sb{x)Sb{x + b/2)Sb{x + b-'/2)Sb{x + g/2) 



(116) 
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